# Hybrid functional calculation for C in the diamond structure
# in a self-consistent approach
# Dataset 1: ground state calculation with WFK output
# Dataset 2-6: calculation of five iterations of HSE06 in the Kohn-Sham basis
# Dataset 7: HSE06 calculation in the planewave basis set
# Dataset 8-13: G0W0 calculations, on top of GGA, HSE06 (Kohn-Sham) and HSE06 (planewaves).
#
ndtset 13
gwpara 2
enunit 1
gw_qprange -14 # Compute correction for all the bands
#gw_qprange -30 # Compute correction for all the bands
prtvol 2
symsigma 0
# Dataset1: usual self-consistent ground-state calculation
# Definition of the k-point grid
ngkpt 2 2 2
nshiftk 4
shiftk 0.0 0.0 0.0 # This grid contains the Gamma point
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
tolvrs 1.0d-15
nband 14
#nband 30
istwfk *1 # Option needed for Gamma
getwfk1 0
# Common to most hybrid/GW calculations
getwfk 1 # Obtain WFK file from dataset 1
ecutwfn 8 # Planewaves to be used to represent the wavefunctions
#ecutwfn 20 # Planewaves to be used to represent the wavefunctions
ecutsigx 8 # Planewaves to be used to represent the exchange operator
#ecutsigx 40 # Planewaves to be used to represent the exchange operator
ecuteps 2
gwcalctyp 25
getqps -1
ixc_sigma -428
pawecutdg 40
gw_icutcoul 6
# Dataset2: Calculation of the 1st HSE06 iteration
optdriver2 4
# Dataset3: Calculation of the 2nd HSE06 iteration
optdriver3 4
getqps3 -1
# Dataset4: Calculation of the 3rd HSE06 iteration
optdriver4 4
getqps4 -1
# Dataset5: Calculation of the 4th HSE06 iteration
optdriver5 4
getqps5 -1
# Dataset6: Calculation of the 5th HSE06 iteration
optdriver6 4
getqps6 -1
# Dataset7: HSE06 stand alone calculation
getwfk7 1
ixc7 -428
# Dataset8: Computation of the dielectric matrix on top of scGGA in the PW basis
optdriver8 3
gwcalctyp8 0
# Dataset9: One-shot G0W0 on top of scGGA in the Kohn-Sham basis
optdriver9 4
gwcalctyp9 0
getscr9 -1
gw_qprange9 0
# Dataset10: Computation of the dielectric matrix on top of scHSE06 in the Kohn-Sham basis
optdriver10 3
gwcalctyp10 20
getqps10 6
# Dataset11: One-shot G0W0 on top of scHSE06 in the Kohn-Sham basis
# (note however that all the bands and k points are considered, in order to perform the correct rotation).
optdriver11 4
gwcalctyp11 20
getqps11 6
getscr11 -1
# Dataset12: Computation of the dielectric matrix on top of scHSE06 in the planewave basis
optdriver12 3
getwfk12 7
gwcalctyp12 20
getqps12 0
# Dataset13: One-shot G0W0 on top of scHSE06 in the planewave basis
optdriver13 4
ixc13 -428 #One should not forget to mention than the starting functional is NOT the one of the pseudopotential ...
getwfk13 7
gwcalctyp13 0
getqps13 0
getscr13 -1
gw_qprange13 0
# Definition of the unit cell: fcc
acell 3*6.7406530878521345 #Same parameters as Shiskin
rprim 0.0 0.5 0.5 #FCC primitive vectors (to be scaled by acell)
0.5 0.0 0.5
0.5 0.5 0.0
# Definition of the atom types
ntypat 1
znucl 6
# Definition of the atoms
natom 2 # There are two atoms
typat 1 1
xred # Reduced coordinate of atoms
0.0 0.0 0.0
0.25 0.25 0.25
# Definition of the planewave basis set (at convergence 16 Rydberg 8 Hartree)
ecut 8 # Maximal kinetic energy cut-off, in Hartree
# ecut 20 # Maximal kinetic energy cut-off, in Hartree
# Definition of the SCF procedure
nstep 250 # Maximal number of SCF cycles
diemac 12.0 # Although this is not mandatory, it is worth to
# precondition the SCF cycle. The model dielectric
# function used as the standard preconditioner
# is described in the "dielng" input variable section.
# Here, we follow the prescription for bulk silicon.
pp_dirpath "$ABI_PSPDIR"
pseudos "C.psp8"
#%%
#%% [setup]
#%% executable = abinit
#%% [files]
#%% files_to_test =
#%% t44.abo, tolnlines=72, tolabs= 6.0e-2, tolrel= 1.1e-0, fld_options = -easy
#%% [paral_info]
#%% max_nprocs = 6
#%% [extra_info]
#%% authors = F. Bruneval and X. Gonze
#%% keywords = GW
#%% description =
#%% Diamond: G0W0 @ scHSE06 calculation. Monitor the direct gap at Gamma.
#%%
#%% First, with the scGW methodology based on a Kohn-Sham basis,
#%% then doing the scHSE06 using the planewave basis, followed by a one-shot G0W0.
#%% The agreement is reasonable with the parameters used in the automatic test,
#%% but can be improved with better parameters (esp. nband), see later.
#%% At the PBE level, the (KS) band gap is 5.231 eV,
#%% At the scHSE06 level, the Kohn-Sham basis delivers 6.834 eV,
#%% while the plane wave basis delivers 6.850 eV.
#%% At the G0W0-scHSE06 level, the Kohn-Sham basis delivers 7.314 eV,
#%% while the plane wave basis delivers 7.344 eV.
#%% The macroscopic dielectric constant (at Gamma) is 9.3698 from PBE,
#%% 5.8309 from scHSE06(KS) and 5.8202 from scHSE06(planewaves).
#%%
#%% These calculations have also been done with better
#%% parameters, in order to observe a better agreement between the KS basis set and the
#%% planewave basis set (ecut 20 ecutsigx 20 nband 30 gw_qprange 30 - note however that ecuteps 2 is low),
#%% at the expense of CPU time..
#%% At the PBE level, the (KS) band gap was 5.661 eV,
#%% At the scHSE06 level, the Kohn-Sham basis delivered 7.340 eV,
#%% while the plane wave basis delivered 7.341 eV.
#%% At the G0W0-scHSE06 level, the Kohn-Sham basis delivered 7.806 eV,
#%% while the plane wave basis delivered 7.807 eV.
#%% The macroscopic dielectric constant (at Gamma) was 7.845 from PBE,
#%% 5.033 from scHSE06(KS) and 5.032 from scHSE06(planewaves).
#%% topics = Hybrids, Susceptibility, SelfEnergy
#%%